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ABSTRACT 

We present a statistical analysis of th e X-ray luminosity o f rotat ion powered pulsars and their 
surrounding nebulae using the sample of Kargalt sev fc Pavlov! (|2008l) and we complement this with 
an analysis of the 7-ray-emission of Fermi detected pulsars. We report a strong trend in the efficiency 
with which spin-down power is converted to X-ray and 7-ray emission with characteristic age: young 
pulsars and their surrounding nebulae are efficient X-ray emitters, whereas in contrast old pulsars are 
efficient 7-ray emitters. We divided the X-ray sample in a young (r c < 1.7 x 10 4 yr) and old sample 
and used linear regression to search for correlations between the logarithm of the X-ray and 7-ray 
luminosities and the logarithms of the periods and period derivatives. The X-ray emission from young 
pulsars and their nebulae are both consistent with Lx oc P 3 /P 6 . For old pulsars and their nebulae 
the X-ray luminosity is consistent with a more or less constant efficiency i] = Lx/E rot ~ 8 x I0~ 5 . 

For the 7-ray luminosity we confirm that L 1 oc \J E rot . 

We discuss these findings in the context of pair production inside pulsar magnetospheres and the 
striped wind model. We suggest that the striped wind model may explain the similarity between 
the X-ray properties of the pulsar wind nebulae and the pulsars themselves, which according to the 
striped wind model may both find their origin outside the light cylinder, in the pulsar wind zone. 
Subject headings: stars: neutron stars - pulsars: general - X-rays: stars 



1. INTRODUCTION 

Despite more than four decades of research, the emis- 
sion from pulsars and their surrounding nebulae is still 
poorly understood. Two important problems are the ori- 
gin of pulsar 7-ray emission and, for pulsar wind nebulae 
(P WNe), the so-called sigma- problem (see the review s 
by IGaensler fc Slandl2006HAjonsl[2r3ol IKirk et al.ll2009h . 
The latter problem derives its name from the sigma- 
parameter, a = B 2 / AirTm c c 2 n, which is the ratio of the 
energy density due to Poynting flux over the particle en- 
ergy density, with B the local magnetic field strength, T 
the Lorentz factor of the pulsar wind, and n the particle 
number density. According to theoretical models, most 
of the pulsar's rotational energy loss rate (or spin-down 
power) E rot is due to Poynting flux, i.e. the pulsar wind 
should have a » 1, but observations of PWNe indicate 
that most of the energy that has been lost is actually con- 
tained by the relativistic electrons/positrons. Somehow, 
the high a flow is converted into kinetic energy some- 
where between the pulsar's magnetosphere and the wind 
termination shock, which converts the radial pulsa r wind 
into an isotropic, relat ivistic particle distribution (jAronsI 
12001 IKirk etaTll2009h . 

A third, and perhaps related, problem is the high wind 
multiplicity factor. The combination of a high magnetic 
and rapid rotation results in a strong electric potential in 
the magnetosphere. This potential will be neutralized by 
charged parti cles that are stripped from the surface of the 
neutron star (fGoldreich fc Julian! I I969t ). The associated 
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charged particle density is 

n GJ = 7 x IQ- 2 B/P cm" 3 , (I) 

with P the pulsar's period and B the local magnetic field. 
A fraction of these particles will escape through open 
field lines, resulting in a particle flux 

N GJ = 2.7 x 10 30 P- 2 Bi2 s~\ (2) 

with B12 the dipole surface magnetic field in un its of 
10 12 G. However, X -ray (e.g. IGaensler et al.l [2002h and 
TeV ()de Ja ger 2007) observations indicate that the num- 
ber of relativistic electrons contained by PWNe turns 
out to be orders of magnitude larger than Nqj, i.e. 
N = kNqj, with the multiplicity factor being n > 500 
for a young pulsar like BI509-58 (|de Jagerl 120071 ) . The 
origin of the additional plasma is likely electron/positron 
pair production in the magnetosphere. The pair produc- 
tion occurs in the presence of the high magnetic fields 
inside the magnetosphere, and requires the presence of 
high energy photons that are either the result of curva- 
ture radiation or inverse Compton scattering. The elec- 
trons that cause the emission are accelerated due to the 
extremely large voltage drop across the open field line s 
(jHibschman fc Aronsll2Ci0ll lHarding fc Muslimovl 12002!) . 
For the inverse Compton scattering seed photons are nec- 
essary that are emitted by the hot polar caps of the pul- 
sar, heated due to the bombardment by accelerated par- 
ticles, or due to the cooling of the young neutron star. 

Despite the many unsolved problems, pulsar research 
has thrived over the last decade thanks to many ad- 
vances in observational techniques and numerical sim- 
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ulations. In particular high energy observations have 
contributed to a wealth of new information on pul- 
sars and PWNe; from high spatial resolution X-ray 
images with Chandra, revealing torii a nd jets (e.g. 
Hester et al.l [200l iHelfand et all l200ll : IBamba et all 



2010blL to a rapid progress in the field of TeV astron- 
omv (jrlinton fc Hofmannll2009t for a review) , which have 
revealed an unexp ectedly large population of very ex - 
tended PWNe (e.g. lde Jagenl2007t iMattana et al.ll2009h . 
This rapid growth in observational data has recently been 
augmented by the GeV 7-ray observatory Fermi, which 
has greatly increased t he sample of 7-ray-emitting pul- 
sars lAbdoiiiI][20T3) . 

Here we present a statistical analysis of two samples of 
rotation po wered pulsars. One is those o f X-ray pulsars 
compiled by [kargaltsc v fc Pavlov! (|2008f ). the other the 
aforementioned sample of -Fermi-detected pulsars. Our 
analysis concentrates on what determines the radiative 
efficiency of pulsars and their PWNe. We report a sur- 
prisingly strong correlation between the X-ray luminos- 
ity of pulsars and their PWNe, which inversely corre- 
lates with characteristic age, at least for young pulsars. 
In contrast, the 7-ray emission correlates po sitively with 
charac teristic pulsar age, as already noted bv lAbdo et al.l 
l(2Ol0h . 

2. A STATISTICAL ANALYSIS 

It is well known that the non-thermal X-ray 
luminosities of pulsars and PWNe are strongly 
correlated wit h the spin-dow n lum in osity E rnt of 
the pulsar dSeward fe Wang] 119881: iVerbunt et al.l 
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19961: IBecker fc Truemperl 119971: IPossenti et al.l I2002t 
Cheng et al.l 120041: IKargaltsev fc Pavlovll2008D . The effi- 



ciency with which the spin-down luminosity is converted 
into X-ray emission is usually indicated by the symbol 



n = Lx/E r o 



(3) 



with rj in the range of rj ~ 10~ 6 — 1. An important ques- 
tion is what determines this efficiency for both 7-ray and 
X-ray emission. For the PWNe rj may provide informa- 
tion on how well the spin-down luminosity is converted 
into relativistic particles, i.e. it is related to the sigma- 
problem. 

For the statistical analysis described here we use the 
X-ray pro perties of pulsars and their nebulae as deter- 
mined by Kargaltscv fc Pavlovl (|2008l ). This sample is 
based on Chandra observations. The Chandra- ACIS in- 
strument that was used has a poor timing resolution, 
so all pulsar luminosities are a combination of pulsed 
and unpulsed emission. However, given the high spa- 
tial resolution of Chandra, the X-ray flux from the pul- 
sar could be accurately separated from the X-ray flux 
from the PWN. The X-ray luminosities were derived from 
the 0.5-8 keV fluxes, corrected for interstellar absorp- 
tion and using the distance estimated listed in the paper. 
Kargaltsc v fc Pavlovl (|2008t ) list only the non-thermal X- 
ray luminosities, so the contribution of thermal X-ray 
emission, if present, was ignored. The uncertainties in 
the distance are the largest source of error in the lumi- 
nosities. In general, the distance estimates may have 
errors of order ~ 2, resulting in luminosity errors of or- 
der 4. Of course distance errors affect the pulsar and 
PWNe luminosities in the same way. A source of er- 



ror for the PWNe luminosities may be underestimated 
flux contributions from l ow surface br i ghtness emission 
at large radii. Indeed, IBamba et al.l (|2010af) recently 
reported the detection of low X-ray surface brightness 
structures around several PWNe by Suzaku. The frac- 
tion of the luminosity in these low surface brightness is, 
however, not more than a factor of 2 in luminosity. So all 
together, the logarithm of the luminosities may contain 
errors of the order of one decade. 

We included in our study a statistical analysis of the 
7-ray luminosities, based on the pulsed 7-ray emission of 
pulsar reported by lAbdo etHTl (|2010( ) . We made sure 
that the 7-ray luminosities of pulsars that were com- 
mon to both the X-ray and 7-ray samples were based 
on the same distance estimate s, namely those adopted 
bv IKargaltsev fc Pavlovl (|2008l ). For both the X-ray and 
7-ray samples we excluded millisecond pulsars, and omit- 
ted pulsars without good distance estimates. 

In figure [1] we show the correlation of the X-ray 
luminosity efficiency 77 of pulsars and PWNe versus 
the characteristic spin-down age r c = P/2P. We 
also included data points for two recently discovered 
PWN e; those surrounding t he high energy pulsar J14003- 
6326 dRenaud et al.1 J2010l) and AXP 1E1547.0- 5408 
(|Vink fc Bambal 120091 1 We did not include these in our 
statistical analysis, although their properties are consis- 
tent with the general trends we report below. 

Figure Q] also shows the 7-ray efficiency, as ob- 
tained from the pu lsed 7-ray luminosities determined by 
lAbdo et al.1 <f2010T> . This fi gure reveals the trend that 
young pulsars appear to have higher X-ray efficiencies 
than old pulsars. Moreover, the behavior is similar for 
the luminosities of the pulsars and the PWNe. 

Figure [2] shows the timing properties of the pulsars in 
the two samples used for our statistical analysis in the 
P — P diagram. It shows that on average the detected 
7-ray pulsars seem to be somewhat older than the X- 
ray detected pulsars/PWNe, consistent with the trend 
that older pulsars are less efficient in emitting X-rays and 
more efficient in emitting 7-rays. Note, however, that an 
interpretation of this diagram is far from straightforward, 
because the detectability of a pulsar also depends on its 
distance and, for X-rays, on the interstellar absorption 
column. 

A trend of decreasing X-ray luminosity w ith r c 
was reported befor e by IBecker fc Truemperl (|1997f> : 
IPossenti et al.l (|2002f) , but figure [1] reveals it to be due 
to a trend in the efficiency, not in the overall spin-down 
power. For pulsars with r c > 2 x 10 4 yr the efficiency 
appears to be more or less constant. In contrast, young 
pulsars are not so e fficient in producin g 7-ray emission, 
as already noted bv I Abdo et al.l (|2010l ). The spread in 
the data points for a given characteristic is of order of 
1-2 decades. As discussed above, the uncertainties in the 
luminosity can explain about 1 decade of this spread, but 
intrinsic variations, and for the 7-ray emission, beaming 
properties are expected to contribute to the spread as 
well. 

In order to investigate the correlations between 77 and 
r c further, one has to avoid rj as an independent variable, 
since both 77 and r c are derived from combinations of the 
pulsar's period P and period derivative P; the spin-down 
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T c . The data has been taken from Kargaltscv & Pavlov (2008) and is based on the X-ray luminosity in the 0.5-8 keV ba nd as measured b y 
Chandra. The open blue squares indicate the 7-ray efficiencies in the fOO MeV to 100 GeV band based on Fermi data (Abdo ct al. 2010). 
The s tar-like symbols indicate two recently discovered pulsars (black) and PWNe (red) which we re not part of the IKargaltsev k Pavlov 
(l2008h sample: the AXP 1E1547.0- 5408 (r ch = 1447 vr lVink k Bamba|[2009l: ICamilo et al.1120071 ). and PSR J14003-6326 {r c h = 23.7 kyr 
IRenaud ct al. 2010). For these two pulsars/PWN the reported unabsorbed fluxes were converted to a 0.5-8 keV band luminosity using the 
distances quoted in the papers. 



luminosity is given by 



E r 



inn = 4tt 2 i 



p3< 



(4) 



with / w 1.4 x 10 45 gem 2 the neutron star's moment of 
inertia. For that reason we base our regression analysis 
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Fig. 2. — The period and period derivative o f the pulsars in the 
X-ray IjKargaltsev fc Pavlovl l200l ) and T-rav ijAbdo et al.l I2010T) 
samples. The symbols match those in Fig. [l] Star-like symbols 
represent AXP 1E1547.0- 5408 (upper right) and PSR J14003-6326. 
The small dots represen t pulsars drawn from the ATNF catalogue 
( Manchester ct al. 20(3). 



on the logarithrrQ of the X-ray/7-ray luminosity, as the 
quantity to be explained, and the independent variables 
log P, log P as the p rinciple input variables for the model 
(c.f. iPossenti et alf[2002l ). 

As explained above the errors are dominated by sys- 
tematic errors, mostly due to difficulties in estimating 
distances. For that reason we used the unweighted least 
square method, which means that the errors in the best- 
fit parameters are based on the variance of the residuals. 
The disadvantage is that we do not have an intrinsic 
goodness of fit statistic. However, we can compare two 
hypotheses using the F-test statistic, defined as: 



F = 



X^(log£» - login) 2 /{m-2 - mi) 
EjOg^ - \ogL 2i ) 2 /{n - m 2 ) 



(5) 



here Li denotes the observed (7-ray /X-ray) luminosity 
and we use Lu to indicate the expected value of Li based 
on the best-fit parameters, and the subscript 1 and 2 to 
indicate two different best-fit models, m is the number 
of degrees of freedom and n is the sample size. The 
probability that the improvement in the sum of squared 
residuals is by pure chance is given by the F(rri2 — mi, n — 
m,2 )-distribution. 

From the above it is clear that we fitted the linear 
relation: 

logL = a + 61ogP + clogP. (6) 
We compare this relation, using the F-test, with a func- 

1 We use log to denote the 10-bascd logarithm. 
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tional form log L — log E rot = a. This is essentially as- 
suming that r\ (equation [3]) is constant. In addition we 
fit the relation 

logL = a + blogE rot , (7) 

which corresponds to L oc E b rot oc P b /P 3b . This 
is a functional form that has been used before (e.g. 
iKargaltsev fc Pavlovl l2C)08f) . and that seems to wor k par- 
ticularly well for the 7- ray emission from pulsars ([Aronsl 
[1991 lAbdo et~aIll20T0h . 

Since figure Q] suggests that young pulsars are more 
efficient X-ray emitters than old pulsars, we divided the 
X-ray sample in an old and a young population. As the 
precise characteristic age that should be chosen is a bit 
arbitrary, but lies somewhere between 10 4 yr and 5 x 
10 4 yr, we decided to divide the sample in two more or 
less equally sized samples. This put the cut at r c = 
1.7 x 10 4 yr. 

The result of the regression analysis is listed in Table [TJ 
whereas figure [3] shows the results of applying the best-fit 
solutions to the different samples. The best-fit relations 
for young pulsars and PWNe indicate that the X-ray 
luminosity does not scale with E rot , since in that case 
we would expect for equation ([6]) b ~ —3 and c ~ 1, 
whereas we find b = —6.2 ± 0.7 and c = 3.3 ± 0.4 and 
b = -5.2 ± 0.7 and c = 2.8 ± 0.4 for PWNe and pulsars, 
respectively. This conclusion is based on both the best- 
fit parameters and their inferred errors, as on the F-test. 
The F-test, in fact, indicates that equation ^ provides 
a better fit than equation ([3]) at the 99.99% confidence 
level. Interestingly, the values for b and c for pulsars and 
PWNe are consistent with each other. 

For older pulsars the F- values indicate that the X-ray 
emission efficiency is just as well described by equation 
([3]) with r] w 8 x 10~ 5 , as by a linear dependency on 
logP and log P. Also this result is true for both the 
pulsar X-ray luminosity as for the PWN luminosity. 

The least square fits for the Permz-sample shows that 

the relation L 7 oc \/ E rot is indeed a good description 
of the data. For equation ([7]) we find b = 0.51 ± 0.10, 
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best-fit models using equation [6] The symbols indicate: black 
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open triangles: luminosity of olds pulsars; blue open squares: 7-ray 
pulsars. 



close to the expected 1/2, and for equation ^ we find 
b = —1.6 ± 0.4 and c = 0.49 ± 0.18, whereas a scaling 

with \J E rot implies b = —3/2 and c = 1/2, in excellent 
agreement. The F-tests confirm this, and indicate that 
a constant value for r] 7 can be rejected at the 99.9996% 
confidence level. 

Finally, one may wonder whether the combined 7-ray 
and X-ray luminosity is more closely correlated with 
E rotl than the 7-ray and X-ray luminosities individu- 
ally. Unfortunately, this is hard to determine from the 
present data sets as there are only 15 pulsars that the two 
samples have in common, with most of them being older 
pulsars. In fact, the combined X-ray/7-ray luminosity of 
all pulsars is dominated by the 7-ray emission, with the 
exception of the Crab pulsar. It is, therefore, not sur- 
prising that the best-fit relation between Lx,^ indicates 

' 1 /2 

Lx +7 oc E rot , just like for 7-ray luminosity. However, 
as the number of 7-ray-detected pulsars will grow in the 
near future, it may be good to investigate the total radia- 
tive output from rotation powered pulsars in more detail, 
especially around r c ss 10 yr, where the radiative output 
changes from X-ray to 7-ray domination. 

3. DISCUSSION 

We presented a statistical study of the X-ray and 7- 
ray properties of rotation powered pulsars, with the aim 
of finding what trends underly the efficiency with which 
spin-down luminosity is converted into high energy radi- 
ation. We started our analysis by showing that young 
pulsars are efficient X-ray emitters, but poor 7-ray emit- 
ters, whereas it is the other way around for old pulsars. 
A trend of low X-ray lumin osity as a function of charac - 
teristic age wa s repo rted by IBecker fc Truemperl (|1997|) ; 
IPossentie t al. (2002), but it is shown here that it is re- 
lated to the efficiency with which spin-down luminosity 
is converted to X-ray emission, and not just due to an 
overall decline in spin-down luminosity. 

In order to be as general as possible in finding trends 
in X-ray/7-rayluminosities our main results are based 
on a regressi on analysis of X-ray lumin osity of the pul- 
sars/PWNe (IKargaltsev fc Pavlovll2008D and 7-ray lumi- 
nosity of pulsars (|Abdo et al.H201Q[) versus the pulse pe- 
riod and its derivative. This is a very generic method, 
which includes the possibility that the luminosity de- 
pends solely on spin-down luminosity. Thus encompass- 
ing other possible dependencies, such as on spin-down 
luminosity. Given the trend noted in figure Q] we divided 
the X-ray sample in two equally sized samples of young 
(t c < 1.7 x 10 4 yr) and old pulsars. 

Our statistical analysis produced two new findings: 1) 
the X-ray luminosity of pulsars and their surrounding 
PWNe appear closely correlated, with Lx :Psr /pwne oc 
j>3 j p6 f or y 0un g pulsars, whereas for old pulsars the best 

fit gives Lx tPsr /pwne oc P I P 35 , the latter being close to, 
and statistically indistinguishable from, a constant X-ray 
emission efficiency 77; 2) young pulsars are more efficient 
X-ray emitters than old pulsars, and have a different de- 
pendency for Lx on P and P. 

We also confir m the findings by I Aronsl (|1996h : 
lAbdo et all (|2010l) that 1) the 7-ray luminosity is well 

described by L 7 oc \/ E rot ; 2) younger pulsars are less ef- 
ficient in producing 7-rays. The latter trend is, therefore, 
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the reverse of the X-ray luminosity (figure [IJ . 

The question is what these findings reveal about pul- 
sar X-ray emission mechanisms. First of all, the connec- 
tion between the X-ray luminosities of pulsars and their 
PWNe may be surprising, given that the luminosities of 
the PWNe are, in general, affected by both the pulsar 
wind properties and the environment of the pulsar. In 
particular for young PWNe the luminosity is possibly 
affected by its interaction with a surrounding supernova 
remnant (SNR). The reverse shock of the SNR will during 
a certain phase of the SNR evolution compress the PWN, 
which naively may be assumed to result in a brightnening 
of the PWN. However, a recent study by Gclfand et al. 
(2009) showed the behavior of the PWN luminosity to 
be more complex. Their study indicates that the com- 
pression by the reverse shock leads to a brightening of 
the radio luminosity, but at the same time to an almost 
total quenching of the X-ray luminosity for a brief period 
(for the specific model they calculated this happened be- 
tween 18-30 kyr). Apart from this brief phase, the X-ray 
luminosity traces the spin-down luminosity surprisingly 
well, with an expected fluctuation of about 0.3 decades 
in rj. 

Given that apparently the X-ray luminosities of the 
pulsars and their PWNe seem to trace each other, and 
have similar dependencies on the period and period 
derivative, one cannot easily use the best-fit functions 
for young pulsars to derive a correlation of X-ray lumi- 
nosity with some well known physical pulsar property, 
such as E rot oc P/P 3 , B p oc VPP, B lc oc pi/2/p5/2 

(the magnetic field at the light cylinder), uqj oc yj P/P, 

or N GJ oc \J P/P 2 oc \J E rot . For young pulsars the 

X-ray emission is poorly fit with a dependency on E h rot 
(Tabic CD). This suggests that the X -ray luminosity 
from pulsars and PWNe may require a more complex 

model than the 7-ray luminosity, with its L 1 oc \J E rot 
scaling. The best-fit formula does, however, explain 
the strong dependency on t c for young pulsars, since 
L x <xP 3 /P 6 <xE rot /(PT?). 

Secondly, the similar behavior of the pulsars and PWN 
X-ray luminosities suggest there is a physical connec- 
tion between the two. It is usually assumed that the 
pulsar X-ray emission originates in the magnetosphere, 
whereas the X-ray emission from the nebulae comes from 
outside the termination shock. These are two distinct 
regions, which are separated by a region that encom- 
passe s the pulsar light cylinder and the so-called wind- 
zone (jKirk et al.ll2009h . the region in which the pulsar 
wind is formed. 

It is tempting to speculate that the connection between 
the X-ray emission from the pulsars and PWNe may have 
something to do with the pair multiplicity. Models of 
pair creation in the magnetosphere indicate that the pair 
multiplicity is a function of the characteristic age (e.g. 
iHibschman fe Aronsl l200ll fi gure 6). This would mean 
that both the X-ray emission from the pulsar and from 
the PWNe are somehow proportional to the multiplic- 
ity. This is not a completely satisfactory explanation, 
because the X-ray synchrotron depends on the total en- 
ergy contained by the pairs, not just by the total number 
of particles. 



Another explanation for the similar behavior of pulsar 
and PWN luminosities is offered by the st riped wind the- 
ory ((CoroniiCIlll EiEHin[lQQllIP09). According to 
this theory the alternating magnetic fields generated by 
obliquely rotating pulsars leads to reconnecting magnetic 
fields in the wind zone. This transforms magnetic energy 
into kinetic energy, thereby changing the p ulsar wind 
from a high to a low a outflow. According to Ki rk et all 
(2002) the zone in which this heating and acceleration 
occurs could be the location of pulsed X-ray emission. 
Such a model makes it easier to explain why there is a 
connection between the pulsar and PWN luminosities, as 
in both cases the energy is generated in the wind zone. 
The higher pair multiplicity of young pulsars may be of 
additional importance as there are simply more particles 
available to be accelerated. 

Finally, there is one other issue to consider, namely, 
why young pulsars do not seem to be efficient in gen- 
erating 7-rays. The fact that it is the other way 
around for the X-ray emission, perhaps indicates that 
the pair multiplicity in young pulsars is so high that 
the electric fields in the m agnetosphere are shorted out 
(jHibschman fc Arons 2001). It may also be worthwhile 
to investigate in more detail what the role of the seed 
photons for inverse Compton upscattering in young pul- 
sars is. The Compton scattered photons (and in some 
cases curvature radiation generated photons) are con- 
verted in pairs, and form th e basis of the pair multiplicity 
(jHibschman fc Aronsll200l"h . The soft seed photons orig- 
inate from the neutron star surface, and may the result 
of polar cap heating due to accelerated electron/positron 
beams, or they are the result of the cooling of the neutron 
star. The latter would induce a dependency on neutron 
star age, which could explain why young pulsars behave 
differently than old pulsars. 

We have only briefly mentioned the TeV detected 
PWNe. The TeV luminosities of PWN e do not seem 
to co rrelate with spin-down luminosity (Mattana et al. 
I2009T) . The probable reason is that the TeV emission 
from PWNe is not so much determined by the current 
energy production of the pulsar, but reflects the time 
integrated energy input from the pulsar wind, which in 
itself scales with the initial spin-period E(t = 0). In X- 
rays the lifetime of the X-ray synchrotron emitting elec- 
trons is short compared to the pulsar lifetime, whereas 
in TeV (and in the radio) this is not the case. 

We can therefore conclude that different parts of the 
electromagnetic spectrum inform us about different as- 
pects of pulsars. The 7-ray emission and pulsed radio 
emission is likely related to what happens in the mag- 
netosphere, whereas TeV and radio emission from the 
PWNe inform us about the time integrated properties 
of pulsars. Based on our statistical study we have added 
here the suggestion that the X-ray emission from the pul- 
sar, like the X-ray emission from the PWNe, may inform 
us about what happens in one of the least understood 
regions surrounding the pulsar: the wind zone. 



JV is supported by a Vidi grant from the Nether- 
lands Science Foundation (NWO). RY is supported by 
grant-in-aid from the Ministry of Education, Culture, 
Sports, Science, and Technology (MEXT) of Japan, 
No. 19047004, No. 21740184, No. 21540259. We 



6 



J. Vink et al. 



have made use of the ATNF online pulsar catalogue: http://www.atnf.csiro.au/research/pulsar/psrcat 



REFERENCES 



Abdo, A. A. et al. 2010, ApJS, 187, 460 
Arons, J. 1996, A&AS, 120, C49 

Arons, J. 2009, in Astrophysics and Space Science Library, Vol. 
357, Astrophysics and Space Science Library, ed. W. Becker, 
373 

Bamba, A., Anada, T., Dotani, T., Mori, K., Yamazaki, R., 

Ebisawa, K., & Vink, J. 2010a, ApJ, 719, L116 
Bamba, A., Mori, K., & Shibata, S. 2010b, ApJ, 709, 507 
Becker, W. & Truemper, J. 1997, A&A, 326, 682 
Camilo, F., Ransom, S. M., Halpern, J. P., & Reynolds, J. 2007, 
ApJ, 666, L93 

Cheng, K. S., Taam, R. E., & Wang, W. 2004, ApJ, 617, 480 
Coroniti, F. V. 1990, ApJ, 349, 538 
de Jager, O. C. 2007, ApJ, 658, 1177 

Gaensler, B. M., Arons, J., Kaspi, V. M., Pivovaroff, M. J., 

Kawai, N., & Tamura, K. 2002, ApJ, 569, 878 
Gaensler, B. M. & Slane, P. O. 2006, ARA&A, 44, 17 
Gelfand, J. D., Slane, P. O., & Zhang, W. 2009, ApJ, 703, 2051 
Goldreich, P. & Julian, W. H. 1969, ApJ, 157, 869 
Harding, A. K. & Muslimov, A. G. 2002, ApJ, 568, 862 
Helfand, D. J., Gotthelf, E. V., & Halpern, J. P. 2001, ApJ, 556, 

380 

Hester, J. J., Mori, K., Burrows, D., Gallagher, J. S., Graham, 
J. R., Halverson, M., Kader, A., Michel, F. C, & Scowcn, P. 
2002, ApJ, 577, L49 



Hibschman, J. A. & Arons, J. 2001, ApJ, 554, 624 
Hinton, J. A. & Hofmann, W. 2009, ARA&A, 47, 523 
Kargaltsev, O. & Pavlov, G. G. 2008, in AIP Conf. Series, Vol. 
983, 40 Years of Pulsars: Millisecond Pulsars, Magnetars and 
More, ed. C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi, 
171-185 

Kirk, J. G., Lyubarsky, Y., & Petri, J. 2009, in Astrophysics and 
Space Science Library, Vol. 357, Astrophysics and Space 
Science Library, ed. W. Becker, 421 

Kirk, J. G., Skjaeraasen, O., & Gallant, Y. A. 2002, A&A, 388, 
L29 

Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, 
AJ, 129, 1993 

Mattana, F., Falanga, M., Gotz, D., Terrier, R., Esposito, P., 
Pellizzoni, A., De Luca, A., Marandon, V., Goldwurm, A., & 
Caraveo, P. A. 2009, ApJ, 694, 12 

Possenti, A., Cerutti, R., Colpi, M., & Mereghetti, S. 2002, A&A, 
387, 993 

Renaud, M., Marandon, V., Gotthelf, E. V., Rodriguez, J., 
Terrier, R., Mattana, F., Lebrun, F., Tomsick, J. A., & 
Manchester, R. N. 2010, ApJ, 716, 663 

Seward, F. D. & Wang, Z.-R. 1988, ApJ, 332, 199 

Verbunt, F. et al. 1996, A&A, 311, L9 

Vink, J. & Bamba, A. 2009, ApJ, 707, L148 



TABLE 1 

Best-fit solutions 



X-ray^ 7" ra >i 
Youngj_ 01c£ 



Fit formula 






Parameter 


PWNe 


PSRs 


PWNe 


PSRs 




(1) log L — log 


;E ro t 


= a 





-3.26 ± 0.33 


-3.61 ±0.23 


-3.94 ±0.18 


-4.22 ±0.15 


-1.22 ±0.15 


(2) log L = a - 


f 6 log 


; E rot 


a 


-37.2 ± 12.2 


-26.4 ± 11.1 


-14.2 ± 7.12 


-1.65 ±6.32 


16.5 ±3.6 








b 


1.91 ±0.33 


1.61 ±0.29 


1.29 ±0.20 


0.93 ± 0.18 


0.51 ±0.10 








f/pL 


7.56/1.2% 


4.22/5.6% 


2.1/17% 


0.17/69% 


24.7/0.004% 


(3) log L = a - 


f 6 log 


; P + c log P 


a 


68.9 ± 5.2 


63.9 ±4.8 


43.1 ±3.9 


39.27 ±3.5 


40.1 ±2.4 








b 


-6.18 ± 0.72 


-5.23 ±0.66 


-4.31 ± 0.78 


-3.18 ±0.70 


-1.56 ±0.39 








c 


3.29 ±0.41 


2.84 ±0.38 


1.11 ±0.28 


0.77 ±0.25 


0.49 ±0.18 








F/ps 


8.60/0.015% 


11.90/0.07% 


1.41/27% 


0.46/64% 


11.88/0.02% 








f/pH 


16.85/0.08% 


15.89/0.1% 


0.78/39% 


0.76/40% 


0.017/90% 



Sample size 19 18 28 



a Based on Kargaltscv & Pavlov (2008), but omitting millisecond pulsars and non-resolved PWNe (#38,40). 

b Based on Abdo ct al. (2013), but omitting millisecond pulsars and pulsars without distance estimates. 

c t c < 1.7 X 10 4 yr 

d r c > 1.7 X 10 4 yr 

C A11 quoted errors are rms errors. 

f F-value and associated probability based on a comparison with formula 1 in this table. 
g Based on a comparison with a fit using formula (1) in this table. 
h Based on a comparison with a fit using formula (2) in this table. 



